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ABSTRACT 



O 

^ [ We present new Very Large Array (VLA) observations toward NGC 253 of 

the recombination line H53a (43 GHz) at an angular resolution of 1''5 x 1''0. 
The free-free emission at 43 GHz is estimated to be ~ 100 mJy, implying a star 
^ ' formation rate of ~ 1.3 Mq yr~^ in the nuclear region of this starburst galaxy. 

I A reanalysis is made for previously reported H92a observations carried out with 

angular resolution of 1''5 x 1''0 and 0''36 x 0''21. Based on the line and continuum 
emission models used for the 1''5 x I'.'O angular resolution observations, the RRLs 
H53a and H92a are tracers of the high-density (~ 10^ cm~'^) and low-density 
(~ 10'^ cm~^) thermally ionized gas components in NGC 253, respectively. The 
velocity fields observed in the H53a and H92a lines (1''5 x 1''0) are consistent. The 
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velocity gradient in the central ~ 18 pc of the NE component, as observed in both 
the H53q; and H92q; lines, is in the opposite direction to the velocity gradient 
determined from the CO observations. The enclosed virial mass, as deduced 
from the H53q; velocity gradient over the NE component, is ~ 5 x 10^ Mq in the 
central ~ 18 pc region. The H92q; line observations at high angular resolution 
(0''36 X 0''21) reveal a larger velocity gradient, along a P.A.~ —45° on the NE 
component, of ~ 110 km s~^ arcsec"^. The dynamical mass estimated using the 
high angular resolution H92q; data (~ 7 x 10^ Mq) supports the existence of an 
accreted massive object in the nuclear region of NGC 253. 

Subject headings: starburst galaxies, radio recombination lines, NGC 253 

1. INTRODUCTION 

NGC 253 is one of the nearest (~ 2.5 Mpc) and brightest starburst galaxies, catalogued 
as an SAB(s)c galaxy (de Vaucouleurs et al. 1976). This galaxy has an inchnation of ~ 79° 
with respect to the hne-of-sight with major axis located at a position angle (P.A.) of 51°, 
also containing a bar-like feature tilted by 18° with respect to the major axis (Pence 1981). 
Observations of NGC 253 have been carried out in the radio (Turner & Ho 1985; Ulvestad & 
Antonucci 1997; Mohan, Anantharamaiah, & Goss 2002; Mohan, Goss, & Anantharamaiah 
2005; Boomsma et al. 2005) infrared (Engelbracht et al. 1998), optical (Forbes et al. 2000; 
Arnaboldi et al. 1995) and X-ray (Weaver et al. 2002) wavelengths. Radio observations, 
which are not affected by dust absorption, are an excellent tool to study the structure and 
kinematics of the nuclear region of NGC 253. Observations in the 21-cm line (Boomsma et 
al. 2005) reveal extra-planar motions of HI that occur at a large scale of up to 12 kpc. High 
angular resolution radio continuum observations (Ulvestad & Antonucci 1997) have revealed 
a number (> 60) of compact sources in the central 300 pc of this galaxy, supporting the 
scenario of a massive star formation episode occurring in the center of NGC 253. Nearly 
half of these compact continuum sources are dominated by thermal radio emission from HII 
regions (Turner & Ho 1985; Antonucci & Ulvestad 1988; Ulvestad & Antonucci 1997). The 
radio continuum and radio recombination line emission, observed at high angular resolution, 
have been modeled using different density components for the ionized gas (Mohan, Goss, & 
Anantharamaiah 2005). The emission models suggest the existence of both low (~ 10^ cm~^) 
and high-density (> 10^ cm~^) ionized gas in the central region of NGC 253. On the other 
hand, the most luminous source (5.79-0.39) is unresolved (< 1 pc) at 22 GHz suggesting the 
existence of an AGN in the center of this galaxy (Ulvestad & Antonucci 1997). Observations 
of broad H2O maser line emission (> 100 km s~^) near this radio continuum source has 
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been invoked as further evidence of the presence of a massive object in NGC 253 (Nakai et 
al. 1995). Mohan, Anantharamaiah, & Goss (2002) modeled the VLA continuum and radio 
recombination hne (RRL) emission for the nuclear region of NGC 253 and favor an AGN as 
the source responsible for the ionization. Observations of hard X-ray emission toward the 
core of NGC 253 were also interpreted as evidence of AGN activity (Weaver et al. 2002). 

The bar-like structure was first observed toward NGC 253 in the near-infrared (NIR), 
covering the inner 150" region of the galaxy (Scoville et al. 1985; Forbes & Depoy 1992). 
The existence of the stellar bar is supported by the observed morphology at optical and 
mid-infrared frequencies (Forbes & Depoy 1992; Pifia ct al. 1992). A counterpart of the 
stellar bar in NGC 253 has been found in CO (Canzian, Mundy, & Scoville 1988), HCN 
(Paglione, Tosaki, & Jackson 1995) and CS (Peng et al. 1996). Observations in the RRL 
H92q; (Anantharamaiah & Goss 1996) at an angular resolution of I'.'S x 1''0 reveal a velocity 
field that is discrepant with the CO, CS and HCN observations. Anantharamaiah & Goss 
(1996) proposed that the kinematics observed in the H92q; could result from a merger of two 
counter-rotating disks. The observed H92a and CO line velocity fields were modeled by Das, 
Anantharamaiah & Yun (2001) using a bar-like potential for NGC 253, which is in reasonable 
agreement with the observed H92« line velocity field. However, this kinematical model can 
only reproduce the velocity field of the CO and CS and does not agree with the H92q; RRL 
observations. Based on the discrepancy of the CO and the ionized gas kinematics. Das, 
Anantharamaiah & Yun (2001) proposed that the accretion of a compact object (~ 10^ Mq) 
about 10^ years ago could account for the velocity field observed in the H92q; RRL. Paglione 
et al. (2004) observed the CO emission at 3" angular resolution for the inner region and 
modeled the kinematics of the molecular gas using a bar potential, concluding that motions 
of the CO gas in the central 150 pc are consistent with a bar potential and report evidence 
of the existence of an inner Lindblad resonance (ILR) . 

Previous interferometric observations of RRLs have been made at low frequencies (e.g. 
~ 8.3 GHz, H92a). VLA observations at 1''5 x 1''0 angular resolution were used by Anan- 
tharamaiah & Goss (1996) to study the kinematics of NGC 253; Mohan et al. (2002, 2005) 
used the VLA observations at 1"5 x 1''0 and 0''3 angular resolutions to determine the phys- 
ical properties of the ionized gas in NGC 253. In this paper we analyse the kinematics of 
the ionized gas in the nuclear region of NGC 253 using high frequency RRL observations 
(~ 43 GHz) and the high angular resolution observations (at 0('3) in the RRL H92q;. Also 
we use the VLA observations in the RRL H53q; and the 43 GHz radio continuum, along 
with previously reported H92q; line and 8.3 GHz radio continuum observations (Mohan et 
al. 2005), in order to estimate the physical properties of the ionized gas. This paper is 
complementary to the results summarized by Mohan et al. (2005). Section 2 presents the 
observations and data reduction, while § 3 presents the results for the H53q! and H92q; RRLs. 
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In subsection 4.1, a model for the emission of the RRLs H53q; and H92q; as well as for the 
radio continuum at 43 and 8.3 GHz is presented. Subsection 4.2 analyzes the kinematics for 
the ionized gas in the center of NGC 253 and § 5 presents the conclusions. 

2. VLA Observations. 
2.1. H53q; line. 

The H53a line {uj-est =43309.4 MHz) was observed in the CnD configuration of the 
VLA on 2003, January 18, 19 and 20. Wc used cycles with integration times of 10 min on 
NGC 253 and 1 min on the phase calibrator J0120-270 (~ 0.7 Jy). Four frequency win- 
dows (LOs) were used to observe the RRL H53q!, centered at 42885.1, 42914.9, 42935.1, and 
42964.9 GHz. For each frequency window, the on-source integration time was ~ 2 hrs, using 
the mode of 15 spectral channels with a channel separation of 3.125 MHz (~ 22 km s~^). 
The data calibration was carried out for each frequency window using the continuum chan- 
nel, consisting of the central 75% of the band. The flux density scales were determined 
from observations of J0137+331 (3C48; 0.54 Jy). The bandpass response of the instrument 
was corrected using observations of J0319+415 (3C84; ~ 7.5 Jy). The parameters of the 
observations are summarized in Table 1. In order to track reliably the phase variations in- 
troduced by the troposphere, the calibration of the data was performed correcting for the 
phases in a first step and subsequently correcting for both amplitude and phase. The line 
data were further calibrated using the solutions obtained by self-calibrating the continuum 
channel of each frequency window. The radio continuum images were obtained by combining 
the continuum channels of each frequency window using the task DBCON from AIPS, and 
the self-calibration method was also applied to this combined data. The H53q; line cubes 
and the 43 GHz continuum image were made using a natural weighting scheme and then 
convolved to obtain a Gaussian beam of 1'.'5 x 1''0 (P.A.= 0°). The combination of the differ- 
ent frequency windows was made following a similar method to that used for the H53a line 
observed toward M82 (Rodriguez-Rico et al. 2004): (1) the line data from each frequency 
window were regridded in frequency using the GIPSY reduction package, (2) before com- 
bining the four LOs into a single line cube, the continuum emission was subtracted for each 
frequency window using the AIPS task IMLIN with a zero order polynomial fit based on 
the line free channels, and (3) the four line cubes (after subtraction of the continuum) were 
combined into a single line cube. The total line bandwidth, after combining all the windows, 
is about 150 MHz (1000 km s~^). The line data cube was Hanning-smoothed using the task 
XSMTH in AIPS to reduce the Gibbs effect and the final velocity resolution is ~ 44 km s~^. 
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2.2. H92q; line. 

We also present previously reported observations of the RRL H92a {urest =8309.4 MHz) 
at X I'.'O, P.A.= 0° (Anantharamaiah et al. 1996) and 0"36 x 0'.'21, P.A.= -3° (Mohan et 
al. 2005, 2002 ) angular resolutions toward NGC 253. The H92a RRL images at 1% x 1"0 
angular resolution were produced by combining observations carried out in the B (August 
31 and Sept 01, 1990), C (May 14 and 23, 1988) and D (July 01 and 19, 1988) configurations 
of the VLA. In order to obtain the same HPFW beam as the H53q; line cube, the H92q; line 
cube was made using these 'B+C+D' combined data applying a natural weighting scheme. 
Anantharamaiah & Goss (1996) have already used these 'B+C+D' combined observations 
to analyze the kinematics of the ionized gas in the central 10" of NGC 253 with an angular 
resolution of 1"8 x 1"0, RA.= 10°. The higher angular resolution (0''36 x 0''21, RA.= -3°) 
observations of the H92q; line toward NGC 253 were made with the VLA in the A array 
(July 9 and 12, 1999) and have been previously reported by Mohan, Anantharamaiah, & 
Goss (2002). Recently, Mohan, Goss, & Anantharamaiah (2005) used these H92q; data 
along with observations in the RRLs H75q; and HlGGo; data to model the RRL and the radio 
continuum emission in order to determine the physical parameters of the ionized gas. In 
this paper we use these high angular resolution observations to study the kinematics of the 
ionized gas in the nuclear 5" region of NGC 253. Because of the different spectral line grid 
of the H92a high angular resolution observations and the H92q; 'B+C+D' data, a combined 
dataset 'A+B+C+D' was not produced. 

For the H92a line observations, the phase calibrator was J0118-216 and the bandpass 
calibrator was J2253+161. A spectral mode with 31 channels was used. The continuum 
images were obtained by averaging the data in the central 75% of the total band. The 
continuum data were processed using standard calibration and self-calibration procedures. 
The calibration and self-calibration used for the continuum data were then applied to the 
line data. All the images were made in the AIPS environment. The line images were 
Hanning-smoothed to reduce the Gibbs effect and the velocity resolution is 56.4 km s~^. 
Further observational details are summarized by Mohan, Anantharamaiah, & Goss (2002) 
and Mohan, Goss, & Anantharamaiah (2005). 

3. RESULTS 

Figure 1 shows the radio continuum emission of NGC 253 at 43 GHz with an angular 
resolution of 1('5 x I'.'O, P.A.= 0° {!" ~ 12 pc). The integrated 43 GHz continuum fiux 
density is 360 ± 20 mJy, obtained by integrating the flux density over the nuclear 30" region 
using the task IRING in AIPS. The radio continuum image at 43 GHz shows two radio 
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continuum components, NE and SW, in addition to extended emission (see Figure 1). The 
continuum peak position of the NE component, a(J2000)= OOMT'^SSf 17 ± OfOl, 5(J2000) = 
— 25°17'17''4 ± 0''1, coincides within 0''2 with the position of the compact source 5.79-39.0 
(Ulvestad & Antonucci 1997). 

Figure 2 shows the H53q; velocity-channel images of NGC 253 with an angular resolution 
of 1('5 X l^'O (P.A.= 0°). The H53q; hne emission is detected toward both the NE and SW 
continuum components above a 3(J level (~ 2 mJy). The ionized gas is observed in the H53q; 
line at heliocentric velocities that range from 17 to 345 km s^^. The velocity-integrated 
H53q; line emission (moment 0) is shown in Figure 3 superposed on the moment of the 
H92q; line. There is good correspondence between the integrated line emission of the RRLs 
H53q; and H92q;. In addition, the peak position of the integrated H53q; line emission is in 
agreement with the peak position of the 43 GHz radio continuum image. In the H53q; line 
images, both the NE and SW components are spatially resolved only along the major axis. 

Figure 4 shows the H53q; hne spectrum integrated over the central 10" region of NGC 253. 
By fitting a Gaussian, the estimated central heliocentric velocity is 210 ± 10 km s~^, the 
FWHM of the line is 230 ± 20 km s"^ and the peak line flux density is 21 ± 2 mJy. The 
resulting flt is shown in Figure 4 along with the residuals to the flt. The central velocity 

is in agreement with previous estimates in the optical (225 ± 5 km s~^; Arnaboldi et al. 
1995) and IR (230 ± 10 km s^^; Prada, Gutierrez & McKeith 1998). The velocity integrated 
H53a line flux density determined from our observations is 0.69 ± 0.09 x 10^^'^ W m^^, in 
agreement with the previous measurement of 0.94 ±0.38 x 10~^° W m~^ derived from single- 
dish observations (Puxley et al. 1997). A Gaussian function was also used to determine the 
characteristics of the spectra obtained by integrating over the NE and SW regions. Table 2 
hsts the results for the total integrated H53q; hne emission proflle, as well as for proflles that 
correspond to the NE and SW components. The values hsted for the H53q; hne are peak 
flux density Sl, FWHM, the heliocentric velocity Vnei and the velocity integrated H53q; hne 
emission. 

Figure 5 shows the velocity held (moment 1) of the H53q; line at an angular resolution 
of 1('5 X I'.'O (P.A.= 0°). Figure 6 shows the H92q; velocity held made using the 'B+G+D' 
data of Anantharamaiah & Goss (1996) at the same angular resolution. The velocity field of 
the ionized gas as observed in the H53a line agrees with observations of the RRL H92q! (see 
section 4.2 for a detailed comparison). In the NE component of NGC 253 the red-shifted gas 
is observed toward the NW and the blue-shifted gas toward the SE. In the region located 
S of the radio continuum peak, there is a blue-shifted component that is more apparent in 
the H53q; line than in the H92q; line. A detailed comparison between the line proflles of the 
H53q; and H92q; in this region shows that the H53q; line is broader than the H92q; line by 
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~ 50 km and there is a relative velocity shift between these two RRLs of ~ 30 km s~^. In 

the elongated SW component, the red-shifted gas is located at the SW and the blue-shifted 
gas is at the NE. The velocity gradient was measured, at 1''5 x 1''0 angular resolution, along 
the major (P.A.= 52°) and nearly along the minor (P.A.= —45°) axis for both RRLs the 
H53q; and H92a;. The H53a velocity gradient along the major axis of NGC 253 (measured 
over the NE component) is 12 ± 3 km s~^ arcsec"^, comparable to the corresponding H92q; 
velocity gradient (~ 11 km s~^ arcsec"^). The velocity gradients measured in the RRLs 
H53q; and H92q; (both at l'/5 x l^'O) along the RA .= -45° are 42 ± 8 km s-^ arcsec"^ and 
24 ± 2 km s~-^ arcsec"^, respectively. 

Figure 7 shows the H92q; velocity field at an angular resolution of 0''36 x 0''21 (P.A.= 
—3°). At this angular resolution, the H92q; hne emission is detected only toward the NE 
component with an angular size of ~ 0'/6 (7 pc). Based on these H92q; data, a larger velocity 
gradient of 110±20 km s~^ arcsec"^ is measured along a P.A.~ —45°. This velocity gradient 
is about a factor of four larger than the velocity gradient estimated using the lower angular 
resolution H92q; observations (1''5 x 1''0). The lower velocity gradient measured in the low 
angular resolution image of H92q; (1''5 x I'.'O) is due to a beam dilution effect. By convolving 
the high angular resolution data of the H92q; with a Gaussian beam of 1''5 x 1''0 angular 
resolution, the velocity gradient measured along the P.A.= —45° is consistent with the lower 
angular resolution H92q; data. 



4. DISCUSSION 

4.1. Models for the radio continuum and recombination line emission 

In Table 3, we summarize the 43 GHz continuum flux density measurement along with 
previous measurements (5 — 300 GHz) for the central 30" region of NGC 253. These flux 
density values have been used to determine the relative contributions from free-free, syn- 
chrotron and dust emission. At frequencies < 50 GHz, the relative contributions of free-free 
and non-thermal emission are dominant compared to the thermal dust emission. However, 
at frequencies > 50 GHz, the dust contribution is more significant. The estimated contri- 
bution of the thermal free-free emission at 43 GHz is ~ 140 mJy, while the non-thermal 
emission accounts for ~ 220 mJy These values were obtained assuming that the thermal 
continuum free-free flux density shows Sfree-free{^) oc Using the observed flux den- 

sity measurements of the radio continuum in the range of 5 — 98 GHz, and following the 
procedure used by Turner & Ho (1983), the spectral index for the non-thermal emission is 
(^synchrotron — —0.73 ± 0.06. After Subtracting the free-free and non-thermal emission from 
the total continuum emission over the 5 — 300 GHz frequency range, we obtain the spectral 
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index value for the dust emission adust — 3.9 ± 0.2. Figure 8 shows the contribution from 
thermal free-free, synchrotron and thermal dust emission along with the total radio contin- 
uum over the frequency range 5 to 300 GHz. The thermal frec-frce radio continuum flux 
density at 43 GHz may be used to estimate the ionization rate {^Lyc) from NGC 253 using 
(Schraml & Mezger 1969; Rodriguez et al. 1980), 

^^-^-9oxio^3r^^-*Mr y-'v ^ y-^fpy 

where Tg is the electron temperature, u is the frequency and D is the distance to NGC 253. 
The estimated ionizing fiux at 43 GHz, assuming Tg = 10"^ K and D= 2.5 Mpc is ~ 7 x 
10^^ s^^. This ^Lyc value is a factor of ~ 5 lower than that estimated by Puxley et al. (1997), 
using the H40q; line fiux density and assuming local thermodynamical equilibrium (LTE). 
However, the RRL H40q; line emission mainly arises from stimulated emission, implying a 
smaller value of the SFR. On the other hand, the star formation rate {^ob) can also be 
estimated using the relation ^Lyc — 5.4 x 10^^ x ^qb (Anantharamaiah et al. 2000), 
obtained assuming a mass range of 1 — 100 Mq in the Miller-Scalo initial mass function 
(IMF). The total SFR derived in the nuclear regions of NGC 253 is thus ~ 1.3 yr'^ 

The emission in the RRLs H53q; and H92q; as well as the radio continuum at 43 and 
8.3 GHz were modeled for the two continuum components (NE and the SW) of NGC 253 
observed at an angular resolution of 1''5 x 1''0 (P.A.= 0°). The line and continuum flux 
densities were measured over regions where both the H92a and H53a arc detected. The 
models consist of a collection of two families of Hll regions, each with different electron 
densities. The electron density ranges that were explored arc n,. = 10^ — 10^ cm^^ and n^ = 
10^ — 10^ cm~^, for the low- and high-density gas components, respectively. The contribution 
from non-thermal synchrotron emission that arise from SNR and the possible AGN and free- 
free emission from HH regions were also considered. The RRL emission has been computed 
considering that the population of the atomic levels deviates from local thermodynamic 
equilibrium (LTE). The non-LTE effects result in both internally and externally stimulated 
line emission. The formalism used to compute the RRL and radio continuum flux densities 
at each frequency follows that used by Rodriguez- Rico et al. (2004). In these models, the 
electron density (ug), temperature (Tg), size (so) and the number of Hll regions (Nhii) are 
free parameters. In order to reduce the number of free parameters, the size of each HH 
region has been assumed to be a function of the electron density and the number of ionizing 
photons emitted by the embedded O star by [so nl^^] — U, where C/ is a function of the 
spectral type of the star (Panagia 1973). In these models, early-type 07 stars have been 
used as the source of the ionizing continuum flux for each HII region. The results obtained 
for Ug and Nhii were constrained by the observed line and continuum flux densities as well 
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as the volume of the total line emission, assuming spherical HII regions. 

The models suggest that the thermally ionized gas in the NE component consists of a 
collection of extended (~ 1 — 4 pc) low-density (~ 10^ — 10^ cm~^) HII regions and compact 
(~ 0.01 — 0.06 pc) high-density (~ 10^ — 10^ cm""^) HII regions. In order to reproduce 
the observations, the mass of ionized gas in the low-density component must be a factor 
of ~ 10^ larger than the mass of ionized gas in the high- density component. The ionized 
gas in the SW region also consists of low- and high-density HII regions, characterized by 
Ug ~ 10^ cm^^ and ng ~ 10^ — 10^ cm^'^. Table 4 lists the physical parameters of these HII 
regions: the electron temperature (Tg), electron density (ng), size, emission measure (EM), 
mass of ionized gas, continuum optical depth (tc) at 8.3 and 43 GHz, departure coefficients 
(b„ and /3„ = 1 — {kTf,/huL)d In bn/dn, where n is the quantum number ) for the H53q; 
and H92ci; RRLs, the contribution of free-free emission and the Lyman continuum photons 
rate. The second and third columns list these parameters for the low and high-density 
HII regions on the NE component. The fourth and fifth columns list the corresponding 
parameters for the SW component. Based on these models, the RRL H92a arises mainly as 
externally stimulated line emission from the extended low-density (~ 10^ cm~^) HII regions, 
with /3„ < -20. 

The thermal free-free contribution from both the low- and high-density component to 
the total observed continuum emission ranges from ~ 30% to ~ 80%. The total mass of 
thermally ionized gas in each component (NE and the SW) is ~ 10^ Mq. The Lyman 
continuum emission rate of 5.9 x 10^^ s~^ obtained from the RRL emission models (see 
Table 4) is consistent with 7 x 10^^ s~^, obtained from the continuum emission models 
(shown in Figure 8). The Lyman continuum emission rates for the NE and SW components 
(fisted in Table 4) were obtained only for regions where H53q; fine emission was detected. 
On the other hand, the value of 7 x 10^^ was estimated by integrating over the central 
30" region, which explains the slightly different results. The results obtained for the low- 
density gas component for the NE component of NGC 253 are in agreement with the results 
obtained by Mohan, Goss, & Anantharamaiah (2005) for the 15 pc region, observed with 
higher angular resolution observations (~ 0''3) and using VLA H166q;, H92q; and H75q; line 
data. Even though the mass contribution of the high density HII regions is only ~ 1% of 
the total HII mass, this feature contributes ~ 40% of the total H53q; fine emission. 

As noted before by Mohan, Goss, & Anantharamaiah (2005), high angular resolution 
observations are essential in the determination of the physical parameters of both low and 
high-density HII regions. In addition to high angular resolution, high frequency observations 
(e.g. H53a) are required to determine the properties of the high-density (~ 10^ cm~^) HII 
regions. Non-LTE effects are important, leading to an enhancement of the line emission by a 
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factor of ~ 2. Thus, the results obtained assuming LTE conditions overestimate the number 
of HIT regions (and also the SFR) that must exist in NGC 253. 

4.2. Kinematics 

The previous H92q; line observations at 1''8 x 1''0 angular resolution (Anantharamaiah 
& Goss 1996) revealed velocity gradients of ~ 11 km s"^ arcsec"^ and ~ 18 km s~^ arcsec"^ 
along the major axis and minor axis, respectively. A qualitative comparison of the velocity 
field observed in the H53q; line and the previously reported H92a line velocity field reveals 
kinematical behaviors that are consistent, i.e. the regions with red-shifted and blue-shifted 
ionized gas coincide (see Figs. 5 and 6). The coincidence of the red- and blue-shifted regions 
in the H53q; and H92a velocity fields implies that both the low- and high-density ionized 
gas components rotate in the same sense. 

In order to compare in detail the kinematics of the ionized gas as observed in the H53q; 
and the H92q; lines, we constructed PV diagrams along the major axis (P.A.= 52°) using 
the task SLICE in GIPSY. These PV diagrams are shown in Figure 9; the white line is 
the resulting fit to the velocity gradient of 11 km s^^ arcsec^^. The same procedure was 
used to obtain PV diagrams along the P.A.= —45°, shown in Figure 10; the white hues are 
the resulting fits to the H92q; velocity gradient (24 km s~^ arcsec"^) and the H53q; velocity 
gradient (42 km s~^ arcsec"^). In Figure 10 we also show the H92q; and H53q; spectra 
superimposed at different offset positions from the central source (5.79-39.0); these spectra 
were normalized based on the peak line flux densities. By inspection of the different spectra 
obtained at the negative offset positions (between —O'.'G and — 0''9 in Figure 10), a relative 
velocity shift (~ 30 ± 5 km s~^) is observed between the peak flux density of the H53q; and 
the H92q; lines. Based on the line emission models, these two RRLs trace different density 
components (see section 4.1). The small velocity shift between these two RRLs on the NE 
component suggests that each density component has slightly different kinematics. 

4-2.1. Gaseous bar structure, outflow or an accreted object 

Observations at IR wavelengths have revealed the existence of a gaseous bar in NGC 253 
(Scoville et al. 1985). In a bar potential the gas follows two types of orbits, xl and x2. The xl 
(bar) orbits are those extended along the major axis of the bar and the x2 (anti-bar) orbits 
are those oriented perpendicular to the bar major axis. In the case of NGG 253, the xl and 
x2 orbits would be oriented on the plane of the sky at P. A. of ~ 70° and ~ 45°, respectively. 
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In the H53q; and H92q; RRL images (at l'/5 x I'.'O angular resolution) the orientation of the 
largest velocity gradient is nearly perpendicular to the orientation of the x2 orbits. Since 
the ionized gas on the NE component rotates in an opposite sense compared to the CO 
(Anantharamaiah & Goss 1996; Das, Anantharamaiah & Yun 2001), a simple bar potential 
does not account for the differences observed between the velocity fields of the RRLs (H92q! 
and H53a) and CO. A secondary bar inside the primary bar may be invoked to explain the 
kinematics observed in the center of NGC 253. However, further observations and modeling 
are required to investigate the existence of this secondary bar. 

Weaver ct al. (2002) proposed the presence of a starburst-driven nuclear outflow col- 
limated by a dusty torus, based on X-ray observations of NGC 253. In this model, the 
thermally ionized gas in the center of NGC 253 should be distributed in both a starburst 
ring and a starburst-driven outflow (Weaver et al. 2002). Observations of the RRL H92q; 
toward the starburst galaxy M82 (Rodriguez- Rico et al. 2004), have proben that RRLs may 
be used to study the ionized gas associated with galactic outflows. The largest velocity 
gradient observed in the RRL H92a is oriented nearly along the minor axis of NGC 253 
(P.A.= —45°). Based on this orientation and assuming the H92q; RRL in the NE component 
traces the ionized gas in the outflow, the receding side of this outflow would be on the NW 
and the approaching side on the SE. If this is the case all the observed ionized gas would be 
tracing the outflow, explaining the different rotation sense between the CO and the ionized 
gas. However, it seems unlikely that all the ionized gas is associated with the outflow. 

Das, Anantharamaiah & Yun (2001) propose that the kinematics of the ionized gas 
traced by the H92q; RRL can be explained if there is an accreted object with mass of 
~ 10^ M0. The CO gas which traces the galactic disk of NGC 253 is moving in an opposite 
sense compared to the ionized gas that may be associated with the compact object. Based on 
the H53q; velocity gradient (~ 42 km s~^ arcsec"^) along the minor axis of the NE component 
(~ 30 pc), the inferred dynamical mass is ~ 5 x 10^ Mq. This mass estimate is consistent with 
that of the accreted object proposed by Das, Anantharamaiah & Yun (2001). The existence 
of a compact object is further supported by the higher angular resolution (0''36 x 0''21) H92q; 
observations (Figure 7), revealing a larger velocity gradient (~ 110 km s~^ arcsec^^, at P.A.~ 
—45°) over the central ~ 0''6 (7 pc). This H92q; velocity gradient of ~ 110 km s~^ arcsec"^ 
imphes a dynamical mass of ~ 7 x 10^ M0, similar to the mass determined from the 1'.'5 x 1''0 
angular resolution observations of the RRLs H53q; and H92q;. The ~ 7 x 10^ Mq dynamical 
mass is based on observations over a region a factor of three times smaller than that observed 
in the 1''5 x 1''0 angular resolution images. 

The estimated dynamical mass (~ 7 x 10^ Mq) for the nuclear region of NGC 253 is 
comparable to that of the compact source at the nucleus of our galaxy (~ 4 x 10^ M©, Chez 
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et al. 2005). This mass estimate of ~ 7 x 10^ Mq for the central region of NGC 253 could 

exist in the form of a large number of stars combined with ionized gas and may also contain 
an AGN. If the Lyman continuum photons rate (~ 7 x 10^^ s^^) is mainly due to 05 stars, 
each emitting ~ 5 x 10"^^ s^^, then there must be ~ 10^ 05 stars in the NE component. 
Using Salpeter's initial mass function and a mass range of 0.1 — 100 Mq, the total mass in 
stars in the NE component is ~ 10^ Mq. The ionized gas could be gravitationally bounded 
by the stars and consequently the black hole mass would be < 10^ Mq. The existence of an 
AGN has been proposed from radio continuum observations (Turner & Ho 1985; Ulvestad 
& Antonucci 1997) and RRL H92q; observations (Mohan, Anantharamaiah, & Goss 2002). 
Radio continuum observations (1.3 to 20 cm), reveal that the strongest radio source 5.79-39.0 
has a brightness temperature > 40, 000 K at 22 GHz and is unresolved (< 1 pc, Ulvestad 
and Antonucci 1997). Broad (> 100 km s^^) H2O maser line emission is observed toward 
the nuclear regions supporting the existence of a massive object in the center of NGC 253 
(Nakai et al. 1995). 

In order to account for the different kinematics observed for the ionized and the molec- 
ular gas, two possible scenarios can be proposed: (1) a dense object that is accreted into 
the nuclear region of NGC 253, (2) the ionized gas is moving in a starburst-driven outflow 
and/or (3) a secondary bar exists within the primary bar, as proposed for other galaxies 
(Priedh & Marinet 1993). The accreted object model is supported by the high angular res- 
olution H92q; observations; the estimated dynamical mass of ~ 7 x 10^ Mq is concentrated 
in a < 7 pc region and the ionized gas traced by the RRLs moves in the opposite direction 
compared to the larger scale CO. In the RRLs H53q; and H92q;, we find no evidence that 
confirms the existence of a secondary bar. An S-shape in the velocity field is characteristic 
of a bar (Anantharamaiah & Goss 1996). Thus, for a secondary bar a second S-shaped 
pattern would be observed in the velocity field which is not appreciated in the H92q; velocity 
structure (see Figure 7). However, this scenario cannot be ruled out and higher angular and 
spectral resolution observations are necessary to discern between these three models. 



5. CONCLUSIONS. 

The H53q; RRL and radio continuum at 43 GHz were observed at high angular resolution 
{1'.'5 X 1''0) towards NGC 253. We have also reanalyzed previous observations of the RRL 
H92q; made at an angular resolutions of 1'.'5 x I'.'O (Anantharamaiah & Goss 1996) and 
0!'36 X 0''21 (Mohan, Anantharamaiah, & Goss 2002). 

Based on the 43 GHz radio continuum fiux density and previous measurements at 
lower and higher frequencies, we have estimated the contribution from free-free emission 
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(~ 100 mJy at 43 GHz). Using this value for the free- free emission, the derived SFR in 
the nuclear region of NGC 253 ~ 1.3 M© yr'^ The RRLs (H53a and H92a) and radio 
continuum (at 43 and 8.3 GHz) emission have been modeled using a collection of Hll re- 
gions. Based on the models, the RRL H53a enables us to trace the compact (~ 0.1 pc) 
high-density (~ 10^ — 10^ cm~^) Hll regions in NGC 253. The total mass of high-density 
ionized gas in the central 18 pc is ~ 10^ Mq. A large velocity gradient (~ 42 km s~^ arcsec"^, 
RA.~ —45°) is observed in the H53q; hne. The orientation and amphtude of the velocity 
gradients (angular resolution ~ 1''2) derived using the H53q; and H92q; lines agree. The 
high angular resolution observations (0''36 x 0''21) of the H92q; hne reveal a larger velocity 
gradient (~ 110 km s~^ arcsec"^); this large velocity gradient implies a dynamical mass 
on the NE component (< 7 pc) of ~ 7 x 10^ M©, supporting the existence of an accreted 
compact object. 

The orientation of the H53q; and H92q; velocity gradients does not agree with the CO 
kinematics. The different kinematics observed over a larger region of the disk of NGC 253 
also suggests the existence of an accreted object. The derived dynamical mass (~ 7x 10® M©) 
can be accounted for by a large stellar density and/or the presence of an AGN. The star 
formation activity in NGC 253 may be the result of a merger process of NGC 253 and this 
proposed massive compact object. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. CR and YG 
acknowledge support from UNAM and CONACyT, Mexico. 
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Table 1. Observing parameters for NGC 253 using the VLA. 



Parameter H53q; RRL (43 GHz) 



Right ascension (J2000) 00M7"33fl8 

Dechnation (J2000) -25°iri7^'0 

Angular resolution continuum l'/5 x I'.'O, P.A.= 0' 

Angular resolution hne l'/5 x V.'O, RA.= 0' 

On-source observing duration (hr) 8 

Bandwidth (MHz) 125 

Number of spectral channels 38 

Center YhcI (km s"^) 200 

Velocity coverage (km s~^) 1000 

Velocity resolution (km s~^) 44 

Amplitude cahbrator J0137 + 331 

Phase calibrator J0120 - 270 

Bandpass cahbrator J0319 + 415 

RMS line noise per channel (mJy/beam) 0.65 

RMS, continuum (mJy/beam) 0.2 



Table 2. H53q; and H92q; line parameters for NGC 253. 



H53a H92a 



Sc43'^ S_L AVfWHM ^Hel Sl ■ AVpwHM ^c&.z'^ Si A.V "VhcI Sl ■ AVfwhm 

Feature (mjy) (mjy) (km s"!) (km s"!) (xlO^^'^Wm-^) (mjy) (mjy) (km s"!) (km s^^) (xlO^^^Wm-^) 



NGC 253 NE 134 ± 5 9.0 ± 0.5 200 ± 10 220 ± 5 0.25 ± 0.01 278 ± 10 4.5 ±0.1 190 ± 5 200 ± 2 2.4 ± 0.2 

NGC 253 SW 24 ±2 2.8 ±0.3 130 ± 20 215 ± 10 0.04 ± 0.01 40 ± 4 0.7 ±0.1 160 ± 15 220 ± 10 0.3 ± 0.03 

NGC 253 Total 360 ± 20 21 ± 2 230 ± 20 210 ±10 0.69 ± 0.09 605 ± 10 9.0 ±0.5 190 ± 10 206 ± 4 4.7 ±0.5 



Continuum emission was obtained by integrating over the region where both RRLs H53a and H92a were detected. 
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Table 3. Observed continuum flux densities for NGC 253^ . 



Frequency 
(GHz) 


Flux density 
(mJy) 


5 


1400 ± 70 ^ 


8.3 


1100 ±25 


15 


640 ± 30 ^ 


43 


360 ± 20 


85 


350 ± 50 ^ 


99 


320 ± 30 ^ 


230 


1700 ± 200 e 


300 


4400 ± 1700 ^ 



^Integrated flux density 
measured over the inner 
~ 30" region. 

•^Turner & Ho (1983). 

^Mohan et al. (2002). 

•^This work. 

^Carlstrom et al. (1990). 
^Peng et al. (1996). 
sKrugel et al. (1990). 
^Chini et al. (1984). 
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Table 4. Results from models using a collection of HII regions for NGC 253. 





NE re 


gion 


sw 


region 


Parameter 


Low density^ 


High density^ 


Low density*^ 


High density*^ 


Te (K) 


7500 ± 2500 


7500 ± 2500 


5000 ± 1000 


5000 ± 1000 


Ue (cm~^) 


1.6 ± 1.4 X 10^ 


5.5 ±4.5 X 10^ 


200 ± 100 


6.5 ±3.5 X 10^ 


Size (pc) 


2.4 ±1.7 


0.04 ±0.03 


4.5 ±1.5 


0.02 ±0.01 


EM (cm-6 pc) 


3 ± 2 X 10^ 


50 ± 45 X 10^ 


1.5± 1 X 10^ 


6 ± 4 X lO'' 


Mhh (Mq) 


3 ± 2 X 10^ 


7±6 


5 ± 3 X lO-'^ 


0.4 ±0.3 


Tc (8.3 GHz) 


6 ±5 X 10-3 


32 ±30 


6±5 X 10-^ 


60 ±40 


rc (43 GHz) 


6 ±5 X 10-5 


1.0 ±0.9 


4±3 X 10-5 


1.9 ± 1.3 


bn (H92«) 


0.881 - 0.973 


0.997- 0.999 


0.819-0.887 


0.998 - 0.999 


bn (H53a) 


0.70 - 0.79 


0.96-0.98 


0.64 - 0.65 


0.97-0.98 


Pn (H92a) 


-55 ± 32 


-1.25 ±0.75 


-70 ± 10 


-0.15 ±0.14 


Pn (H53o;) 


-27± 13 


-10 ±3 


-13 ±2 


-5±1 


S//_43^(mJy) 


70 ±30 


19 ±2 


NLyc (XIO^^ g-1) 


5±2 


0.9 ±0.1 



^Both components were used to fit the continuum emission at 8.3 and 43 GHz as well 
as the flux densities in the RRLs H92ci; and H53q;. 

^Free-free continuum flux density at 43 GHz obtained summing the contributions from 
the low and high-density components. 
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Fig. 1. — Radio continuum image of NGC 253 at 43 GHz obtained using the VLA. Contour 
levels are drawn at -3, 3, 6, 12, 24, 48, 96, 144, 192, 240, 288, 336, 384, and 432 times 
the rms of 0.2 mJy beam~^. The cross shows the position of the compact source 5.79-39.0 
(Ulvestad & Antonucci 1997). The angular resolution is 1'.'5 x V.'O, P.A.= 0°. 
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Fig. 2. — Channel images of the H53a line emission toward NGC 253 obtained using the 
VLA. Contours are —3, 3, 4, 5, 6, 7, 8, 9 and 10 times 0.65 mJy beam~^, the rms noise. 
The cross shows the position of the compact source 5.79-39.0 (Ulvestad & Antonucci 1997). 
The synthesized beam (1('5 x I'.'O FWHM, P.A.= 0°) is shown in the first panel. The central 
heliocentric velocity is given above for each image. 
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Fig. 3. — The velocity integrated (moment 0) of the H53q; hne emission (contours) super- 
posed on the moment of the H92q; hne emission (gray scale) toward NGC 253. Contour 
levels are drawn at 5, 10, 15, 20, 30, 90% the peak hne emission of 1.16 Jy beam~^ km s~^. 
The gray scale (moment of H92q;) covers the range 0.18 — 18 Jy beam~^ km s~^. The cross 
(+) shows the position of the compact source 5.79-39.0 (Ulvestad & Antonucci 1997). The 
HPFW, for both the H53a and the H92a hne images is 1''5 x I'.'O, P.A.= 0°. 
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Fig. 4. — Total integrated H53« spectrum from NGC 253, obtained by integrating over the 
central 10". The thick line shows the data, the dashed line shows the Gaussian fit and the 
dash-dotted line shows the residuals to this fit. 
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Fig. 5. — Velocity field observed in the H53q; line from NGC 253. Contour levels are the 
heliocentric velocities at 180, 190, 200, 210, 220, 230, 240, 250, 260, and 270 km 8"^ The 
gray scale ranges from 170 to 280 km s~^. The cross (+) shows the position of the compact 
source 5.79-39.0 (Ulvestad & Antonucci 1997). The HPFW is 1^'5 x l^'O, P.A.= 0°. 
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Fig. 6.— The H92q; velocity field of NGC 253 (Anantharamaiah & Goss 1996), is shown 
in contours superposed on the gray scale image of the same H92q; velocity field. Contour 
levels are drawn for the hehocentric velocities of the ionized gas from 160 to 260 km in 
steps of 10 km s~^ and the gray scale covers the same heliocentric velocity range. The cross 
(+) shows the position of the compact source 5.79-39.0 (Ulvestad & Antonucci 1997). The 
HPFW is 1('5 X l^'O, P.A.= 0°. 
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Fig. 7. — High angular resolution (0''36 x 0''21, P.A.= —3°) heliocentric velocity field (mo- 
ment 1) image of the H92a line emission toward NGC 253 obtained using the VLA A array 
data. Contour levels are drawn at 160, 170, 260 km s~^ The cross (+) shows the position 
of the compact source 5.79-39.0 (Ulvestad & Antonucci 1997). 
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Fig. 8. — Integrated radio continuum emission measured toward the nuclear 30" region of 
NGC 253 at 5, 15 GHz (Turner & Ho 1983), 43 GHz (this work), 85 GHz (Carlstrom et 
aL 1990), 98 GHz (Peng et al. 1996), 230 GHz (Krugel et aL 1990) and 300 GHz (Chini 
et al. 1984). The contributions to the continuum flux density are: free-free {S free- free oc 
ly"'^'^), synchrotron {S synchrotron oc i/~°-^) and dust {Sdust oc i/^"^) emission which are shown 
as indicated. 
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Fig. 9. — Position-velocity comparisons along the major axis (P.A.= 52°) of NGC 253. The 
H53q; RRL PV diagram (color scale) was constructed from the data with angular resolution 
of 1''5 X 1''0, P.A.= 0° and velocity resolution of 44 km s~^. The color scale covers the 
range 0.64 — 6.4 mJy beam~^. The H92q; RRL PV diagram (contours) was constructed 
from the data of Anantharamaiah & Goss (1996) with angular resolution of 1''5 x I'.'O, 
P.A.= 0° and velocity resolution 56 km s^^ Contours are 10%, 20%, 90% of the 
peak intensity 3.9 mJy/beam. The white dashed line shows the fitted velocity gradient of 
11 km s~^ arcsec"^. 





Fig. 10.— Position-velocity diagrams along a P.A.= -45° in NGC 253. The H53a RRL PV 

diagram (color scale) was constructed from the data with angular resolution of 1''5 x I'.'O, 
P.A.= 0° and velocity resolution of 44 km s~^. The color scale is from 0.64 to 6.4 mJy/beam. 
The H92q; RRL PV diagram (contours) was constructed from the data of Anantharamaiah 
& Goss (1996), with angular resolution of 1''5 x 1''0, P.A.= 0° and velocity resolution of 
56 km s~^. Contours are 10%, 20%, . . . , 90% of the peak intensity 3.9 mJy/beam. The white 
dashed hnes marked with 'a' and 'b' show the fitted velocity gradients of 42 km arcsec"^ 
and 25 km s~^ arcsec"^, respectively. The panels show the H53q; (blue) and H92q; (red) hne 
profiles for the central 1''8. Each panel corresponds to a slice taken along the velocity axis 
of this PV diagram e.g. the panel at 0.0 is a shce along the velocity axis at offset position 

o:'o. 



